Spin transition rates in nanowire superlattices: Rashba spin-orbit coupling effects 
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We theoretically investigate the influence of Rashba spin-orbit coupling in a parabolic nanowire 
modulated by longitudinal periodic potential. The modulation potential can be obtained from 
realistically grown supperlattices (SLs). Our study shows that the Rashba spin-orbit interaction 
induces the level crossing point in the parabolic nanowire SLs. We estimate large anticrossing width 
(approximately 117 fieV) between singlet-triplet states. We study the phonon and electromagnetic 
field mediated spin transition rates in the parabolic nanowire SLs. We report that the phonon 
mediated spin transition rate is several order of magnitude larger than the electromagnetic field 
mediated spin transition rate. Based on the Feynman disentangling technique, we find the exact 
spin transition probability. For the case k = 0, we report that the transition probability can be tuned 
in the form of resonance at fixed time interval. For the general case (k ^ 0), we solve the Riccati 
equation and find that the arbitrary values of k induces the damping in the transition probability. 
At large value of Rashba spin-orbit coupling coefficients for (k ^ 0), spin transition probability 
freezes. 



I. INTRODUCTION 

Low dimensional semiconductor nanostructures such 
as quantum dots, quantum wells and quantum wires 
can be formed in the plane of 2 Dimensional Electron 
Gas (2DEG) with the application of externally applied 
gate potentials, have attracted significant interest for 
building robust spintronics logic devices and other ap- 
plications Single electron spins in these nanostruc- 
tures can be manipulated by several parameters such as 
the gate controlled electric fields in the lateral direction 
and externally applied magnetic fields. The Rashba and 
Dresselhaus spin-orbit couplings provide another efficient 
way to control the single electron spins in these nanos- 
tructuresjii^ The Rashba spin-orbit coupling arises due 
to structural inversion asymmetry in the crystal lattice 
along the growth direction^ The Dresselhaus spin-orbit 
coupling arises due to bulk inversion asymmetry in the 
system^ 1 - 

Accurate estimation of the spin transition rate, medi- 
ated by phonons and electromagnetic fields, are of great 
interest for the design of optoelectronic devicesi 7 ' 12 ' 13 
Long spin relaxation rates, approximately 0.85 ms in 
GaAs quantum dots and 20 ms in InGaAs quantum dots, 
have been measured by utilizing several different experi- 
mental techniques such as pulsed relaxation and optical 
orientation methodsi 14 ' 15 In these experiments, it is con- 
firmed that the transition rate is dominated by the spin- 
orbit coupling with respect to the environment & 16 ' 1 7 Be- 
cause of the spin-orbit coupling, the electron spin qubits 
in the nanowire quantum dots localized in a transmis- 
sion line resonator can be manipulated, stored and read 
out with the application of the gate controlled elec- 
tric fields i 12 ! 18 In this paper, we theoretically investigate 
the energy spectrum of the parabolic nanowire modu- 
lated by realistically grown SLs in the longitudinal di- 



rection>i2r~— We focus on the study of the crossing of 
the energy spectrum of the nanowire SLs accounting for 
the Rashba spin-orbit coupling. The crossing point can 
be achieved with the accessible values of the strength 
of the Rashba spin-orbit coupling. We investigate the 
spin transition rate in parabolic nanowire SLs with the 
Rashba spin-orbit coupling under the influence of electro- 
magnetic field radiation , 22 ' 23 phonons^ and Dyakonov- 
Perel (DP)2i mechanisms. We can write the momentum 
as a classical variable in the Dyakonov-Perel mechanism 
under the Markovian process 2 ^ - — and estimate the tran- 
sition probability by utilizing the Feynman disentangling 
technique method ^ 27 ' 28 The DP mechanism corresponds 
to the spin splitting of the conduction band in zinc blende 
semiconductors at finite wave vectors which is equivalent 
to the presence of an effective magnetic field that causes 
the precession of an electron spinj 24 i 25 

The paper is organized as follows: in section ILT1 we de- 
velop a theoretical model that allows us to find the cross- 
ing of the energy spectrum of the parabolic nanowire SLs 
with Rashba spin-orbit coupling. In this section, we also 
develop the theoretical model that allows us to find the 
spin transition rate. In section IIII1 we plot the disper- 
sion relation (see Figs.[T| of the parabolic nanowire mod- 
ulated by longitudinal periodic potential. In Fig. [2] and 
[31 we plot the spin transition rate of the nanowire SLs 
with Rashba spin-orbit coupling via electromagnetic field 
radiation, phonons and Dyakonov-Perel mechanisms. Fi- 
nally, in section ITVl we summarize our results. 



II. THEORETICAL MODEL 

The total Hamiltonian of the quasi one dimensional 
parabolic nanowire formed in the plane of 2DEG with 
Rashba spin-orbit coupling in presence of longitudinal 
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modulation potential can be written a; 



,4,29-32 



H,. 



e„(k), 



H — H xy + H z + Hi 

2^ + 2 m ^ V 

H z =p 2 z /2m + V(z), 

H R = — (Pydx -PxVy) , 



(1) 

(2) 

(3) 
(4) 



where H z is the Hamiltonian of the electrons along z- 
direction and Hr is the spin-orbit Hamiltonian asso- 
ciated to the Rashba coupling. The full descriptions 
of H z and Hr will be discussed shortly. p x and p y 
are the momentum operator, m is the effective mass, 
a is the strength of the Rashba spin-orbit coupling and 
ujq = h/ mt\ is the confinement frequency of the parabolic 
potential with l§ being the oscillator strength, (i = 
x,y,z) are the Pauli spin matrices. e v (k) provides the 
periodic longitudinal modulation potential along x direc- 
tion in the form of^ ' 33 i 34 



e„(fc) = ^(l 



cos kl) 



(5) 



where Ax is the first miniband width and I is the super- 
lattice periodi 20 ' 35 The second term in (fTJ) represents the 
Hamiltonian of the electron along z-direction where V z is 
the asymmetric triangular quantum well confining poten- 
tial along z-direction. Usually, the asymmetric triangu- 
lar quantum well potential can be found by solving the 
Schrodinger-Poisson equations self-consistentlyi 3 ^ The 
potential along z-direction can be chosen as V z = eEz 
for z > and V z — oo for z < 0£1 The ground state 
wavefunction (^q z (z)) of H z can be written in the form 
of Airy function (Ai) a&L ' 36 ' 37 



^oz(z) = lA261q^ 2 Ai (qz + <pi), 



(6) 



where (pi — —2.3381 is the first zero of the Airy function 
and 



q = 



2meE 



1/3 



(7) 



Since [H xy ,H z ] — 0, we use the momentum along z- 
direction as a good quantum number i.e., < p 2 >= 
0.78 (hq) and < z >= 1.56/g to estimate the thickness 
of the 2DEG. The structural inversion asymmetry in V z 
leads to the Rashba spin-orbit coupling (see Eq. @| and 
the Rashba coefficient a can be written as 3 -^ 

h 2 A {2E g 



7#e < E > 
a = ,7fl 



A) 

2mE g {E g + A) {3E g + 2A) 



,(8) 



where A stands for the spin-orbit splitting in the valence 
band and E g is the band gap. For InAs material, we 
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adopt the value 7^ = 110 A . 

Since, [p x ,H xy ] — 0, we consider p x is the good quan- 
tum number and the eigenvalue of p x can be written 



as-;- 



29,30,34 



m de v (k) ml 



dk 



2h 



Ai sin kl. 



(9) 



The Hamiltonians ([2]) and (HJ can be written in terms 
of annihilation and creation operators as 

H xy = (a i a+ - J hojQ + e v (k) (10) 



ami .... ia 
Hr = 5-A1 sin kl <j v H = 

2h 2 v £ V2 



a a x 



11) 



To find the energy spectrum of the above Hamil- 
tonian, it is convenient to rotate the Hamiltonian 
H xy R = exp(—ma x /4:)(H xy + H R )exp(iircr x /4;) so that 
the eigenvalue of p x couples to <j z . The new Hamiltonian 
can be written as H xy R — Hq + Hi, where Hq is the 
diagonal part and Hi is the nondiagonal partj^ 

H = (a^a + ij huj - ^-^sinkl a z + e v (k), (12) 



Hi 



)V8 



a) cr+ + H.c, (13) 



where £ = Ti 2 /(am) is the spin precession length, a± — 
a x ±i<r y and H.c. stands the Hermitian conjugate. From 
Eq. (TlB")) . it is clear that the non-diagonal part couples 
the state differed by one quantum number. In a case, 
where Hq » Hi, we use non-diagonal Hamiltonian Hi 
as a perturbation. Based on the second order perturba- 
tion theory, the energy spectrum of the nanowire can be 
written as 



-n,+l/2 



- I nujQ — - sin kl - 

n 



- e v (k) 



lAi sinkl 



hujo 



iAi sinkl 



,(14) 



-n,-l/2 



huj a + 1 sin kl + e u (k) 



a 



TlLOn 



|Ai sin kl 



TlLOn 



|Ai sin kl 



.(15) 



We now turn to the calculation of spin-flip transition rate 
in parabolic nanowire modulated by periodic potential. 



Electromagnetic field mediated spin transition 
rate 



We apply the formalism of time dependent pertur- 
bation theory to the interactions of parabolic nanowire 
modulated by longitudinal periodic potential with the 
classical electromagnetic radiation fieldi 21 ' 39 The total 
Hamiltonian of the nanowire SLs with electromagnetic 
field radiation is given b y 22 i 40 



H 



H A = A 

m 



H +Hi 
P 



ae 

~^r~A y (7 x . 



(16) 
(17) 



Eq. (I17|) is treated as an external perturbation dependent 
on the position r and time t. The vector potential A (r,t) 
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FIG. 1. (Color online) Dispersion relation: energy cliff, vs k. (a) spin splitting energy increases and becomes maximum at 
k — tt/21. (b) It can be seen that the crossing between the states |0, — 1/2 > and jl,+l/2 > takes place approximately at 
k — 0.19/nm and k = 0.43/nm. Here we find the anticrossing width is approximately 117 ueV. We chose the material constants 
for InAs material as m = 0.0239, huo — 1.3 meV, Ai = 16 meV and I = 5.06 nm. 




FIG. 2. (Color online) (a) Electromagnetic field mediated spin transition rate vs k in InAs nanowire SLs. (b) Phonon mediated 
spin transition rate vs k. The dips in the transition rate (see solid lines) can be seen due to level crossing. Here we chose 
hw = 1.3 meV, eh 14 = 0.54 x 10~ 5 erg/cm, s; = 4.2 x 10 5 cm/s, s t = 2.35 x 10 5 cm/s, p = 5.6670 g/cm 3 , e r = 14.6, 
Ai = 16 meV and I = 5.06 nm. 



of the electromagnetic field radiation is written as 

A(r, t) = J —^e qX b q , x e^ r -^ + H .c, (18) 

where uj q — c|q|, b q \ annihilate photons with wave 
vector q, c is the velocity of light, V is the volume 



of the nanowire, e r is the dielectric constant of the 
nanowire. The polarization e q A with A = 1,2 are 
chosen as two perpendicular induced photon modes in 
the nanowire. The polarization directions of the in- 
duced photon are e q \ = (sin0, — cos0, 0) and e q 2 = 
(cos 9 cos 0, cos 9 sin 0, — sin 0) because we express q = 



4 




FIG. 3. (Color online) Transition probability, 101/2,-1/2 vs time in the InAs nanowire SLs. Here we chose Ai = 72 meV, 
I — 53 nra, huo = 1.3 meV, £0 = 1.5ao and ao = 0.053 nm is the atomic Bohr radius. 



q (sin 9 cos 4>, sin 9 sin 0, cos 9) . The above polarization 
vectors holds the relation as e q \ — e q 2 x q, e q 2 = qx e q \ 
and q = e q i x e q 2. Based on the Fermi Golden Rule, 
the electromagnetic field mediated spin transition rate 
(i.e., the transition probability per unit time) in the 
nanowire modulated by longitudinal periodic potential 
is given by>22 

7f = ~ 72Z / l M 9,A| 2 ^(^ q - £o,-l/2 + £0, + l/2^ 



A=l,2 



(19) 

where the matrix element M q .\ =< n, l/2\HA\n, —1/2 > 
has been found perturbatively. The spin transition rate 
(i.e., (n, -1/2) (n, +1/2)) is given by 



1 a 2 e 2 Ai^sinfcZ 
T\ 47r?i 4 ^e r eoc 3 



1 - 



4£ 2 



(20) 



Here, p is the crystal mass density, V is the volume of 
the nanowire SLs, 6 qce creates an acoustic phonon with 
wave vector q and polarization e Q , where a = I, t\, t% are 
chosen as one longitudinal and two transverse phonon 
modes. Also, Aq a = q\qke,fiijk^ a is the amplitude 
of the electric field created by phonon strain, where 
q = q/g and efiijk = eft.14 for i 7^ k,i 7^ j,j 7^ k. 
The polarization directions of the induced phonon are 
ei = (sin 9 cos 4>, sin 9 sin 0, cos 9), e tl — (sin 4>, — cos </>, 0) 
and et 2 = (cos 9 cos 0, cos 6* sin</>, — sin 0). Based on the 
Fermi Golden Rule, the phonon induced spin transition 
rate in the nanowire SLs is given b y 37 ' 41 
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1 



(eM 2 « 2 (^ 



iki 



70tt?i 



(huo) 2 - ( Al '| infc ' )' 



1 



4 1 

34 



a- 



Ml 



o (huj a ) 2 - ( 



A±l sin kl 



(22) 



B. Phonon mediated spin transition rate 

We now turn to the calculation of the phonon mediated 
spin relaxation rate in parabolic nanowire SLs. Following 
Ref.0, the interaction between electron and piezo-phonon 
can be written a o 17 ' 41 " — 



<h(^) = \ hr^r—e^-^eA cla bt ia + H.c. (21) 



C. Spin- flip transition probability: Feynman 
disentangling technique 

We used the Feynman disentangling technique^ ' 27 ' 28 to 
find the transition probability of electron spins of the 
Hamiltonian associated to the Rashba spin-orbit cou- 
pling. Under the Dyakonov-Perel mechanism ) 24 ' 25 we 
consider the momentum p(t) — mr as a classical variable 
whose dependence on time t is generated by a Markovian 
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process^ In this case, the Rashba spin-orbit Hamilto- 
nian Hr = exp{—iira x / A)H R exp{iira x / A) can be written 
as 



H R (t) 



^fiuiQ cosujot (s + + s_) — sinfc^ s 2 , (23) 



where £ = ?iu;o^o/a is the spin-precession length, £o 
is the spin-orbital radius and s± = s x ± is y . The 
spin operators obey the 5/7(2) algebra, [s + ,s_] = 2s z 



and [s 2 ,s±] = ±s±. The spin evolution operator, 

U(0,t) — T exp(^—i/h J dtH R (t)^j can be exactly found 

by utilizing Feynman disentangling method 2 -^ U(0, t) = 
exp(a(t)s+)exp(b(t)s z )exp(c(t)s-), where a(t),b(t) and 
c(t) are time dependent functions that can be found ex- 
actly. In the disentangled form, the evolution operator 
can be written as 



U(t) 



a(t)s + >T ft'-o [(t^" cos "ot — x) s' + — ^Ai sin kl s'z + ^-fi-uo cosu t s'_] dt' 



(24) 



where 



III. RESULTS AND DISCUSSIONS 



a(t) = -jJ X (t')dt', (25) 
Jx(.t')dt' s ^^ s+ J x (t')dt'_ ^ 



By differentiating Eq. (|21)1) with respect to 'a' and utiliz- 
ing the initial condition s^(0) = s M , we find the relations 
= s + , s' z = s + a + s z , s'_ = s_ — 2s 2 a — a 2 s + . By 
substituting these relations in Eq. (|24|) and equating the 
coefficient of s+ = 0, we find the differential equation in 
the form of 



da 
~di 



•j^ftwo coswoi (l — a 2 ) — a—A\ sinfcZ | 



(27) 

In an analogous way, we can disentangle sq and s_ by 
differentiating Eq. ([26| with respect to 'b' and 'c' respec- 
tively. However, a single function a(t) needed to find the 
transition probability^ The differential Eq. (|27|) can be 
solved exactly for the case k — 0. Thus its solution can 
be written as 



*(*) = 



exp ^— 2i^ sinc^ot) — 1 



exp y— 2i^jr smujQt \ + 



1 



(28) 



The transition probability between the opposite spin 
states can be written as 

w+i/2 ,-i/2 = 1 jjTj^p = sin 2 ^ sinw tj • (29) 

From Eq. [29l it is clear that the spin flip transition prob- 
ability is enhanced with the Rashba spin-orbit coupling. 
At large value of Rashba spin-orbit coupling constant, 
the splitting of the peak value in the transition probabil- 
ity can be achieved due to the fact that the periodicity 
of the propagating waves in the crystal lattice changes 
with a. For the general case k ^ 0, we solve the Riccati 
Eq. (|27[) numerically to find the transition probability. 



In Fig. [lja), we have plotted energy vs k of InAs 
parabolic nanowire SLs in presence of the Rashba spin- 
orbit coupling. We find the Kramer's type degeneracy 
point at k = ±mr/l with n = 0, 1, 2 • • •. For the arbitrary 
values of k, the degeneracy can be lifted and we find the 
spin splitting energy. The splitted energy difference in- 
creases with the Rashba spin-orbit coupling strength and 
becomes maximum at k = tt/21. Note that this maxima 
point is also maxima of p x (see Eq. [9]). In Fig. Q] (b), we 
investigate the level crossing point associated with the 
spin states \n, —1/2 > and \n + 1, 1/2 >. The crossing 
takes place at k 0.19, 0.43/nm. Here the width of the 
anticrossing is approximately estimated as 117fieV. 

We now turn to the key results of the paper: the 
spin transition rate via electromagnetic field radiation, 
phonons and Dyakonov-Perel mechanismi 17 i 23 ~ — 

In Fig. [5J we plot the spin transition rate vs k between 
the spin states |0, — 1/2 > and |0, +1/2 > under the in- 
fluence of electromagnetic field radiation and phonons. 
It can be seen that the phonon mediated spin transition 
rate is several order of magnitude larger than the value 
mediated by electromagnetic field radiation because the 
electromagnetic field mediated transition rate vanishes 
like sin kl whereas, the phonon mediated spin transition 
rate vanishes like (sinfcZ) 5 (see Eqs.l2"0"landl2"2j). The dips 
in the spin transition rate (see Fig. [21 solid lines) can be 
seen due to level crossing. The dips in the spin transi- 
tion rate can be tuned with the application of the Rashba 
spin-orbit coupling. The cusp like structure in the spin 
flip rate was reported by the authors in Refs. and 
|45| for quantum dots system. For the case of parabolic 
nanowires modulated by longitudinal periodic potential, 
the dips found in the spin transition rate is new. 

Finally, in Fig. [31 we plotted the transition probability 
vs time. In Fig. [21(a) and (b), we show that at k = 
(see solid lines), the spin transition probability can be 
tuned with fixed time interval in the form of resonance. 
It can also be seen that the spin flip transition probabil- 
ity can be enhanced with the Rashba spin-orbit coupling. 
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We study the influence of the inclusion of the wavenum- 
ber (k). We find that there is a superposition effect be- 
tween p x and p y in the evolution of the spin dynamics 
(see Eq. [M] and |2"T)) which induces the damping effect in 
the transition probability with respect to time. At large 
value of a, the third term of Eq. ff2T|) dominates over the 
first and second term and we find the spin freezing in the 
transition probability (see dashed line in Fig. [3Jb)). Re- 
cently, similar type of results have been shown in Ref. |32| . 

IV. CONCLUSION 

In Figs.Q] we have demonstrated that the level cross- 
ing in parabolic nanowire SLs can be achieved with the 
accessible values of the strength of the Rashba spin-orbit 
coupling. The crossing point can be tuned to the lower 
values of k and vice versa with the application of the 
Rashba spin-orbit coupling. In Fig. [3J we have shown 



that the phonon mediated spin transition rate is several 
order of magnitude larger than the electromagnetic field 
mediated spin transition rate. The dips in the spin transi- 
tion rate can be found due to level crossing. Based on the 
Feynman disentangling technique method, in Fig. [3l we 
have shown that the transition probability can be tuned 
in the form of resonance at k = 0. For the arbitrary 
values of k, we have shown that there is a superposition 
effect which induces the damping in the transition prob- 
ability. At sufficiently large values of the Rashba spin- 
orbit coupling coefficients, the spin transition probability 
freezes for the arbitrary values of k. 
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